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ABSTRACT 

The Newtonian n-Body Problem is modified assuming positive inertial masses 
but different sign for the interacting force which is assumed with the pos- 
sibility of two different signs for the gravitational masses, according to the 
prescription two masses with same sign attract one to the other, two masses 
of different sign repel one to the other. As in electrostatics the signed mass 
is called charge. The inertial mass is always positive. The two body problem 
behaves as the similar Coulomb problem of charged particles with two equal 
charges. The solution is a central configuration with almost same behavior 
that the Newton two-body problem for hyperbolic orbits. The 3-Body prob- 
lem was found with collinear solutions. The four body case of charged central 
configurations has only the planar [TJ and collinear solutions. 
Keywords Few-Body problem. Physical central configurations. Charges of 
both signs 
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1 Introduction 



In this paper we present the study of central configurations of few particles (2, 
3, 4) that obey Newton's three laws of motion and Newton's gravitational 
force law [2], [3], assuming the possibility of gravitational charges of both 
signs, according to the prescription (which is opposite to that in electrostatics 
for electric charges where G would have a negative value) that the force 
between two masses with charge of the same sign is attractive, while the 
force between two masses with charge of opposite sign is repulsive. More 
explicitly, we assume the equations of motion 

d 2 Ti v-^ Geien 



' "E^fe-*,). V, (1) 



where ij denotes the position vector of particle j in M 3 , rrij is its positive 
mass, G is the positive constant of universal gravitation, ry = |rj — r ^ | is 
the distance between particles j and /, and e,- is the charge of particle j such 
that rrtj = \ej\, with two possible choices of sign for the charge Cj. 

In the following we quote some fundamental equations of classical me- 
chanics that are essential in Physics [2], [3]. 

Equation of motion (1) expresses Newton's second law equating the posi- 
tive inertial mass times the acceleration to Newton's gravitational force. This 
force obeys the action-reaction law or Newton's third law: the force vector 
that particle j exerts on particle I is of equal magnitude and opposite sign 
to the force that particle I exerts on particle j. As a consequence, the sum 
over all j of the various equations of motion is the null vector 

Vj 

The few bodies' masses mi, m2,- are positive and generally different in 
value, but some may be equal. 
The total mass is 

m = J2 m i' ( 3 ) 

Vj 

The center of mass position is defined as 

c = ^S m ^- ( 4 ) 

Vj 
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Equation d2J) implies 

d 2 c . s 

which asserts that the center of mass moves with constant velocity. 
With no loss of generality we assume in this paper that 

^ = 0, c = 0, J2 m ^ = - ( 6 ) 

The center of mass is thus at the origin of the system of coordinates rj. 
From equation ([1]) the conservation of total energy E follows, namely 

1 v — \ dvj drj v — ■> GciGi ,_, ,_, 
-> m,— r---r--> — — = E. (7) 

2 ^ 3 dt dt ^ ru V ; 

In this expression the first term on the left hand side is the kinetic energy, 
involving the positive inertial masses, while the second term is the potential 
energy, which depends on the gravitational charges. 

Equation ([!]) also implies conservation of angular momentum 



4 V rrijTj x = , (8) 
dt^ J 3 dt w 

Vj 

which again contains the inertial masses. 



2 The integrable two body problem 

In this section the positions of the two particles are written in terms of the 
relative position r = r 2 — ri as 

m 2 m 1 
ri = r , r 2 = — r . (9) 

m m 

The differential equations of motion are 

£ = v, ^ = Gm±r. (10) 
dt dt | r- 1 3 v 1 

With the constants of motion of energy and angular momentum written in 
terms of, the specific energy 

E=-\v\ 2 + Gm^- (11) 
2 1 1 r v ' 



3 



and the areal velocity 

g = rxv = c/(0,0,l), (12) 

where g is twice the magnitude of the areal velocity. This last implies the 
orbit (r and v) is in a plane orthogonal to the constant vector along g. Polar 
coordinates: r, ip, in this plane give us 

r = r (cos x/j, smip, 0) , v = r(cosij;, sin^, 0) + rip(— sin^, cosi[), 0) , (13) 

where the dot on a letter denotes the time derivative, and 



g = r 2 ip, E = -(r 2 + r 2 ip 2 + Gm-) = -(r 2 + + Gm-) . (11) 
2 r 1 r r 



Dividing the second equation in (10) by the first equation in (14) one has 

dw 1 dv Gm . 

(cosip, smip, 0) , (15) 



dip ip dt g 
which is integrated into 

Gm 



fj 



sin?/', — cos?/', 0) + h , (16) 



where h is a constant vector of integration, the Hamilton vector [3], in the 
plane of the orbit. Vector v — h traces a circle in velocity space with center 
at h and radius — . 

It is useful define the constant of motion called the Laplace-Runge-Lenz 
vector defined here in terms of other constants of motion as 

1 -h x g = -Lv x g + - = e(l, 0, 0) , (17) 



Gm Gm r 

which defines the direction of one coordinate axis in the plane of the orbit. 
Projecting vector r in this direction (in polar coordinates) lead to the orbit 
equation 

er cos?/' = h r , (18) 

Gm 

2 

which is a hyperbola with e > 1, the eccentricity and with p = the latus 
rectum. The specific energy becomes E = 
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3 Central configurations of few charged masses. 
Are there equilateral few body central con- 
figurations? 

We define a central configuration as one in which the particle positions, 
defined up to rotation around the center of mass and uniform dilatations, 
obey the equation 

s "* r i = £ " r i) ' V ^ (19) 

¥i lj 

where B is the same quantity for all the particles. In the right side we have 
the gravitational force acting on particle j due to the other charges. Using 
the fact that the sum of forces is the zero vector, we recover the third element 
of hypothesis fl6]): the origin of coordinates r,j is at the center of mass. 

For the particular case of the previous section with just two particles 
the allowed hyperbolic motion is at any moment in a central configuration 
with B = ^Tp. In the following, one considers the cases with three and four 
particles. 

The force of the right hand side of f ll9p is the gradient of the potential 
energy 

E^^H^E — • (20) 

The left hand side of (fT9l) is related with the gradient of the total moment 
of inertia, I, 

d_l 

3 Vfc 

which may be written in terms of the relative distances as 

1 = \^2 m k r k-r k = ^^2m k mir 2 kl . (22) 

Vfc k^l 

It follows that the condition for central configurations of few particles in 
three dimensions may be expressed in terms of derivatives with respect to 
the relative distances rh in the form 

= amirrij , (23) 



mjTj = £ mkrk ' Vk ' ( 21 ) 
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where a is a quantity containing B and G. The left hand side of this equation 
is the derivative of the potential energy with respect to the square of the 
distance rfj, while the right hand side is proportional to the derivative of the 
total moment of inertia with respect to the same variable rfj. 

Theorem 1 

Equilateral central configurations do not exist for three or four particles with 
gravitational charges of different sign (the proof for four particles is at [1]). 

Proof: Equation (23) can not possibly be satisfied for all combinations 
of indexes because the left hand side takes different signs and the right hand 
side has the sign of a because the m's are positive. One concludes that there 
is no equilateral central configuration for particles with different charges. 
Thus, the Lagrange equilateral triangle solution [5], and the Lehman-Fillies 
equilateral tetrahedron solution [6] valid for positive charges are not a central 
configuration for the case of positive and negative charges. □ 

4 Collinear Three-Body and planar Four-Body 
central configuratios for masses of different 
charge 

Consider now the question of possible collinear/planar configurations for 
Three/Four-Bodies. The modified Dziobek equations for collinear/planar 
central configurations of positive and negative charges are 

= amirrij + XSiSj , (24) 

where A is a new parameter and Sk are the directed distances/areas of the 
segments/triangles having at their vertexes the two/three particles differ- 
ent from k. With respect to equation ( l23l) . the additional term takes into 
account the zero area/volume restriction necessary for collinear/planar con- 
figurations: SiSk is a constant times the derivative with respect to rf k of 
the square of the area/volume of the triangle/tetrahedron formed by the 
three/four particles, which is expressed as the square of the Heron's equa- 
tion/as the so-called Cayley-Menger determinant. A similar equation was 
obtained by Dziobek [7], but although the two terms on the right hand side 
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of this equation are the same as in Dziobek's paper, in the left hand side the 
charges are replaced by masses. 

Note the invariance of the modified Dziobek's equations with respect to 
a sign change of all the charges, as well as their invariance with respect to a 
sign change of the directed distances/areas 

e j ^ ~~ e j i > —Sj . 

In order to prove that these equations are equivalent to the equations 
defining the central configurations, we need the collinear/planar conditions 

5>,- = 0, 5>r, = 0. (25) 

Vj Vj 

Substituting the left hand side of (1241) by the right hand side in ffl9l) . one 
obtains 

4 

BrrijVj = G^^(ri — Vj)[amimj + XSiSj] = —GammjYj , (26) 

¥i 

where we used the conditions (EJ) that the center of mass is at the origin, and 
the planar configuration properties (j2"5|) . Thus from ( }2"6"|) we obtain that 

B = -Gam . (27) 

In addition, one has the property 

J2 sis k ri = ^211 s * s s ^ - 2r ' • r i + r ?) = > ( 28 ) 

1=1 3=1 1=1 3=1 

where we use the collinear/planar configuration conditions (I25p . Note that 
the terms with j = I in these summations are zero; therefore if we multiply 
both sides of equation ff24|) by rf- and we sum over all I and j with I ^ j one 
obtains 

7^ = a m i m i r l ■ ( 29 ) 

In this case a has not a definite sign because of the presence of charges with 
both signs. 
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5 Three-Body collinear central configurations 

We ordered the coordinates of three particles as x\ < x 2 < x 3 .The Dziobek 
equations (|24|) for the collinear Three-Body Problem central configurations 
are 

eie 2 1 



mim 2 (x 2 — xi) 3 

e 2 e 3 1 
m 2 m 3 (x 3 - x 2 ) 3 

e 3 ei 1 



where is 



Theorem 2 



m 3 m\ (x 3 — X1) 3 

(x 3 - x 2 )(x 2 - xi)(xx - x 3 ) 



-L. R 1713 


(30) 


x 2 - X\ 


a + (3 


(31) 


x 3 - x 2 


± a 7712 


(32) 


X\ - x 3 



(3 = X- 



mim 2 m 3 



(33) 



It is impossible to have a collinear Three-Body central configuration of dif- 
ferent charges with the middle charge equal to one of the other charges. 

Proof: Assume the different charge is e 3 = — m 3 . The previous equations 
fl30ll32l) are in such a case 



1 


{X 2 






-1 


(x 3 


- a; 2 ) 3 




-1 


(x 3 


- Xx) 3 



a + 

a + f3 
cr + /3 



m 3 



x 2 - Xi 
m\ 

x 3 - x 2 
m 2 



Xi - x 3 

Canceling a between (I34p and (136 p we have 



+ 



(x 2 - Xi) 3 (x 3 - Xi) 3 







m 2 



m 3 



+ 

x 3 — x\ x 2 — Xi 



(34) 
(35) 
(36) 

(37) 



that implies f3 > 0. 

Canceling a between ( )36|) and (|35|) we have 



(x 3 - x 2 ) 3 (x 3 - Xi) 3 



m 2 



mi 



+ 

X 3 — X\ x 3 — x 2 



(38) 



that implies (3 < 0, which is a contradiction. □ 

On the contrary, when the middle charge is the one that has the opposite 
sign, the previous equations (I30H32P are 

: a + Z?^ 3 — (39) 
x 2 - Xi 

- a + f3^^- (40) 

- a + Z?^ 2 -. (41) 
Xi - x 3 

Canceling a and /3 among these three equations we have the equation for 
collinear central configurations 

m 2 (x 3 - x 1 ) 2 (x 2 - xi) 3 + m 3 (x 2 - x 1 ) 2 (x 3 - x 2 ) 3 + m 3 (x 2 - a;i) 2 (x 3 - xi) 3 = 

m 2 (x 3 -a;i) 2 (a;3-a;2) 3 + mi(a;3-a;2) 2 (a;3-a;i) 3 +^i(a;3-a;2) 2 (a;2-si) 3 • (42) 
Particular collinear central configurations are 

x 2 - Xi 





1 


{x 2 


-xi) 3 




1 


(x 3 


- x 2 ) 3 




1 


(x 3 


- Xi) 3 



x 3 - x 2 
x 2 — X\ 3 



2 , mi = 34 m 2 = 10 , m 3 = 5 . 
mi = 105 m 2 = 84 , m 3 = 20 . 



x 3 - x 2 2 

6 Four-Body planar solutions 

Write equation f[2~4"j) divided by the product of masses mirrij as 

= (r + XAiAj , (43) 

mirrij r\- 

where Aj = Sj/rrij denotes, as in [8J, weighted area. Noting that the square 
of e ' gj equals positive one, this equation may also be written as 



1 \ 3 

eiej 1 
mirrij r tj 



o + XAiAj , (44) 
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A similar equation is the basic tool to compute central configurations in refer- 
ence [5] from the given weighted areas Aj. The difference is the replacement 
of distance ry by a sort of charged distance 

rnirrij , . 

rij — y J -rij. (45) 

e iej 

The algorithm presented in j8] is also useful for the present case if allowance 
is made for these positively or negatively charged distances in the numerical 
computation. In that paper the algorithm to compute planar central con- 
figurations starting from the weighted areas was applied to several examples 
with positive charges. Since the areas of the triangles having its vertices at 
the positions of the particles were computed by Heron's formula in terms 
of the square of this distances, they are not modified by the extra factor 
{miirij) /{eitj). At the end of the present paper the modified algorithm is 
used to obtain some numerical examples of central configurations combining 
it with the new algorithm recently proposed to compute central configura- 
tions from given masses [9] based on a new system of coordinates [10] which 
will be described in the next section. We refer the reader to those references 
for a more detailed account. 

Substracting term by term equation (|4U]) with subscripts l,j and I, k, and 
with subscripts n, j and n, k, with all subscripts different, yields 



e z e 7 - 1 eie k 1 



mirrij rfj mim k rf k 



XAi(Aj - A k ) , (46) 



- = ^A n (Aj - A k ) . (47) 
m « m i r nj ™ n mk K k 

Elimination of A between these two equations gives the fundamental relation 

( ejm k e k mj\ ( ejm k e fc mA 
S n ei ^ = Sie n ^ . (48) 

V r lj r lk J \ r nj r nk J 

The same equation was obtained directly from two of the defining equations 
for planar central configurations (1T§|) . A similar calculation may be found in 
reference [8]. 
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Theorem 3 



A planar central configuration of four particles with one charge opposite in 
sign to the other three leads to a concave configuration with the different 
sign charge in the convex hull of the other three (the same proof is at pQ). 

Proof: Assume that the charge of particle 1 is of opposite sign to the 
charge of the other three particles. Using equation (|48p with j — 1, we note 
that the quantity in the two parenthesis have the same sign and therefore 

f >0, f >0, f >0. (49) 
03 04 02 

From the above equation we prove with no loss of generality that the sign 
of the charge of any particle may be made to coincide with the sign of the 
corresponding area. 

S 2 >0, S 3 >0, S 4 >0, 5i<0.D (50) 

Theorem 4 

A planar central configuration of four particles with two pairs of particles of 
different sign leads to a convex configuration with charges of the same sign 
located at the ends of the two diagonals (the same proof is at []]). 

Proof: Consider now the case in which particles 1 and 2 have charges 
of one sign and particles 3 and 4 of the opposite sign. Again using equation 
( 14~8]) with charge j of one sign and charge k of the opposite sign, the quantity 
in the parenthesis have the same sign, so that 

|<0, !<0, §U0, f <0. (51) 
03 04 03 04 

With no loss of generality we take 

S 1 >0, S 2 >0, S 3 <0, S 4 <0.D (52) 

As a consequence, the quotient S^/e^ may be considered to be always 
positive. 
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Theorem 5 



To have a Four-Body central configuration determined by Dziobek-like equa- 
tion ( 124"]) . parameter A is positive (the same proof is at [1]). 

Proof: Divide both members of equation (1241) by the product of charges 
titj to yield 

-i = + A^ . (53) 

Choosing two of these equations, with subscripts say l,j and I, k, such that 
the term containing a in each is of opposite sign, adding member by member 
we obtain 

Since the quantities in the parentheses are both positive, we have proved that 

A > . □ (55) 



7 New coordinates in the Four-Body Prob- 
lem 

This section reviews the main ideas and results of the new four-body coordi- 
nates of [10J, slightly expanded at a few spots but condensed to the minimum 
necessary. This coordinate system has been used to compute several central 
configurations with gravitational charges of both signs. 

We transform from the inertial referential system to the frame of principal 
axes of inertia by means of a three dimensional rotation G parameterized 
by three independent coordinates. In addition to this rotation, three more 
coordinates are introduced, as scale factors R±, R 2 , R3, which are three 
directed distances closely related to the three principal moments of inertia 
through 

h = M^2 + Rl) , h = n(Rl + R{) , and I 3 = n{R\ + R 2 2 ) , (56) 
where \x is the mass 

/ mi m 2 m 3 m^ 
ym 1 + m 2 + m 3 + 777.4 
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The first rotation changes nothing and after the scale change the result- 
ing four-body configuration has a moment of inertia tensor with the three 
principal moments of inertia equal. The second rotation G' does not change 
this property. 

The cartesian coordinates of the four particles, with the origin at the 
center of gravity, written in terms of the new coordinates are 



X\ x 2 
Vi 2/2 
z\ z 2 





(5? 



where G and G' are two rotation matrices, each a function of three indepen- 
dent coordinates such as the Euler angles, and where the column elements 
of the constant matrix 



d\ CL 2 <^3 Q"l 

E = I bi b 2 b 3 6 4 

C\ c 2 c 3 c 4 



(59) 



are the coordinates of the four vertexes of a rigid orthocentric tetrahedron 
[TTj . having its center of mass at the origin of coordinates, namely: 



airrii + a 2 m 2 + a 3 m 3 + a 4 m 4 = , 
bimi + b 2 m 2 + 6 3 m 3 + 6 4 m 4 = , 
qmi + c 2 m 2 + c 3 m 3 + c 4 m 4 = , 



(60) 



We introduce the mass matrix 

/ mi 



V 



M 





m 2 








m 3 




\ 




m 4 / 



(61) 



An equivalent condition for having three equal moments of inertia for the 
rigid tetrahedron is expressed as 



EME 




(62) 



The system of coordinates for measuring the G' rotation can be chosen in 
various ways, from which we prefer to use the same coordinates as in reference 
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[PJ, namely, particle with mass mi along coordinate axis 3, the other three 
in a plane parallel to the coordinate plane containing axes 1 and 2 but that 
does not include the particle of mass mi; the particle with mass m 2 on an 
orthogonal coordinate plane that contains the first particle and the center 
of mass, and the other two particles on a line that is parallel to coordinate 
axis 1 and perpendicular to the coordinate plane containing the first two 
particles. Particle 1 thus has coordinates 



mim 



Particle 2 has coordinates 




K6 2 ,c 2 )=IO,W " (m3+m i - t /^^l. (64) 



Particle 3 has coordinates 



rri2(m — mi) ' y (m — mi)m 
(a 3 ,b 3 ,c 3 ) 



111714 / 111712 / 



m 3 (m 3 + m4) , y (m 3 + m 4 )(m — mi) ' y (m — mi)m 
Particle 4 has coordinates 

(04, 64, C4) = 



(65) 



/im3 / /xm2 / //mi 



(66) 



^4^3 + 1714) ' y (m 3 + m^m — mi) ' V{m — mi)m j 
Note that 63 = 64 and 02 = 03 = C4, as they should. 

This rigid tetrahedron is the generalization of the rigid triangle of the 
Three-Body problem with the center of mass at the orthocenter discussed 
previously in [12]. The same triangle was used with different purposes by C. 
Simo [T3] . 

Our coordinates are now adapted to the important and old subject [7] of 
planar configurations, with the four particles in a constant plane. Since the 
z-component of the four particles equals zero, the first rotation is just by one 
angle in the plane of motion and the scale factor associated with the third 
coordinate is zero, namely 





x 2 


%2, 


X,i 


Vx 




y?> 
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cos if> — sin if) 
sin ip cos if> 
1 

This equation simplifies to 



Rt 
R 2 




G 



CL\ Q>2 ^3 ^4 

bi b 2 b 3 6 4 

Ci C2 C3 C4 



(67) 



X2 %3 £4. 
V\ V2 V3 Va 



cos if) — sin if> 
sin -0 cos if} 




(68) 



CME T G' 




(69) 



in terms of six degrees of freedom. 

The corresponding expression for the four directed areas in terms of these 
coordinates is 

/^\ 

S 2 
S3 

where C is a constant with units of area divided by mass. Note that Equa- 
tions (1251) are satisfied from this expression of the directed areas since sub- 
stitution of equations (1581 and (1591 in equations (1251 . and application of 
equation (1621) . yields an identity, independent of coordinates. 

For G' two rotation angles are needed, for which we chose those required 
to express the unit vector in spherical coordinates 



G' 



sin 6 cos (f> 
sin 6 sin 
cos 9 



(70) 



where 9 and are the spherical coordinates determining this vector. Given 
the four masses, the four directed areas are functions, up to a multiplicative 
constant C depending on the choice of physical units, of this unit vector 
direction only. 

The non-collinear planar central configurations are characterized in our 
coordinates by constant values of the G' matrix and of the constant value of 
the ratio R1/R2, which are not arbitrary, but they are determined by three 
independent quantities as discussed in the following. 
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2 




Figure 1: Stereographic projection of the hemisphere of the two angles motion 
of the orthocentric tetrahedron. The great circles represent the positions 
where three particles are collinear. The four spherical triangles are concave 
open sets labeled by the particle at the interior of the triangle. The spherical 
quadrilateral open sets correspond to convex configurations with the same 
order that the neighboring triangles. The isolated points are at the angles 
where a charged central configuration has been computed in pp. The values 
of the masses are mi = 10, m2 = 13, = 15, WI4 = 17. The sign of the 
charges is opposite for the position inside a concave spherical triangle set 
and apposite for the two particles along each diagonal when the particle is 
inside a rectangular convex region. 
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From ( 1tj§j) follows that in a planar solution the weighted directed areas 
are expressed as 



A 3 
\A 4 J 



CE T G' 




sin 9 cos < 
CE T I sin sine 
cos 9 



(71) 



Therefore, the weighted directed areas are up to a normalization factor equal 
to the third rotated coordinate of the rigid tetrahedron. In terms of the 
vectors (I63H66I) and the angles 9 and cf> this equation is expressed simply as 



Aj = C (a,j sin 9 cos <fi + bj sin 9 sin <p + Cj cos 9) . 

Choosing C = \J (m — m±) / /J, we have explicitly 

m — ml 



(72) 



A 1 



y/mim 



cos # . 



A, 



A, 



A 9 = - 



mi 

— cos 9 




ni 



mi 

— cos 9 




m 



mi m 3 + m 4 . . 

— cos 9 + \ sin 9 sin < 

m V m 2 

sin 9 sin + 
sin # sin </> - 



m 2 



m 3 + m 4 



mAm — mi) 
sin 9 cos ( 

m 3 m 3 + m 4 ) 



m 2 



^3 + m 4 



mAm — mi) 
sin 9 cos < 

m^(m 3 + mA) 



(73) 

Note that a sign change of the unit vector (47) produces a simultaneous 
sign change in these four A,'s which does not to give a different central 
configuration. Therefore, it suffices to consider only the hemisphere < 9 < 
vr/2. 



8 Computing central configurations for posi- 
tive and negative charges 

In the case of positive charges, since the lengths and masses are defined 
up to arbitrary units, with no loss of generality we assumed in [8] that the 
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parameter a equals plus one. However, central configurations with charges 
of different sign obtained numerically always yield a negative a, so that for 
those cases we used a = — 1 as follows 

Ji*-r]* = -l + \A j A k (j^k). (74) 
rrij rrik J 

We recall that in the paper by Pina and Lonngi [8] the assumption that 
the directed weighted areas are known as four given constants was made. 
The equation which corresponds to (51) then gives the distances as functions 
of the unknown parameter A. Through them, the areas of the four triangles 
become functions of A, that should satisfy restrictions fl25|) for a planar solu- 
tion. These restrictions allow in many cases to determine the value of A and 
hence the values of the six distances and the four masses. This is an implicit 
way to deduce planar central configurations with four masses. 

In contrast, in reference [T], as well as in reference [9], we assume that the 
four masses are known from the beginning. The four weighted areas are then 
determined by expressions ( !73|) in terms of the two tuning variables 9 and 4>. 
Particular values of these two angles determine the four constants Aj, up to 
a multiplicative factor, which in turn produce a central configuration with 
computed distances and masses. The computed masses are in general not 
equal (or proportional) to the starting values used to build the orthocentric 
tetrahedron. The two angles are then tuned until a numerical match is pro- 
duced between the given and the computed masses. The distances between 
particles, computed for this central configuration, correspond to the given 
masses. 

For the arbitrarily chosen values of the masses m x = 10, m 2 = 13, m 3 = 
15,m4 — 17 seven central configurations were found in [1] which are repre- 
sented in Figure 1. 



9 Collinear Four-Body central configurations 

The collinear case is very similar to the planar case. Assume the particles 
are ordered in the line with coordinates: r\ < r 2 < r 3 < r 4 .In this case the 
configuration is determined by three distances between particles a = r 2 — r\, 
b = r 3 — r 2 , c = r 4 — r 3 . Using the fact that the coordinates are defined with 
respect to the center of mass: m\T\ + 777.2^ + 7Ti 3 r 3 + 777,47-4 = 0, it is possible 
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r(7) > r(13 

r(6) > r(13) 

Figure 2: The set of the four-body collinear configurations in the central circle 
by the stereographic projection of the hemisphere of the two angles position 
of the orthocentric tetrahedron. The great circles represent the positions 
where two particles have the same coordinate. Each point on the hemisphere 
represent a different collinear configuration. The particular permutation is 
determined by the inequalities of the coordinates associated to two masses 
on both sides of the great circles. Each of the interior of the twelve spherical 
triangles represent the set of different collinear configurations with the same 
permutation order. The values of the masses are m\ = 20, m 2 = 13, = 
7, m 4 = 6. The inequalities between coordinates in figure are labeled with 
the numerical value of the corresponding masses. To simplify compilation 
the r's in figure are not italic fonts. 
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to write the positions r~ in terms of the masses and the distances a, b, c 







( 


—a — (a + 


b) -(a + b + c) \ 


/ 


mi 
m 


\ 


r 2 






a 


-b 


-(6 + c) 




m 2 
m 










a + b b 





— c 




mz 
m 




\u J 




[a 


+ b + c b + 


c c 


o / 


V 


nij 
m 


/ 
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(75) 
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is skew-symmetric with zero determinant. It is expressed in terms of the 
orthogonal vectors 
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In order to have a collinear central configuration vector with entries rj, or the 

F ■ 

parallel F with entries will be perpendicular to two linearly independent 
vectors orthogonal to the vectors in this equation. Two particular, linearly 
independent vectors are any two of the vectors 
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which are linearly dependent according to 

feki + ck 3 = (a + b + c)k 2 



(78) 
(79) 



Theorem 6 

It is impossible to have a collinear Four-Body central configuration with a 
different charge outside the other three charges. 
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Proof: Assume the different charge is at r 4 , then the forces in terms of 
masses and distances are 

F\ m 2 m 3 m 4 



m x G a 2 (a + b) 2 (a + b + c) 2 
F 2 mi m 3 m 4 

m 2 G a 2 b 2 (b + c) 2 

F 3 mi m 2 m 4 

~m 3 G = "(a + 6) 2 ~~ 1? " I 2 



(81) 
(82) 



Fa mi m 2 m 3 



m 4 G (a + b + c) 2 (b + c) 2 c 

Interior product of F with the vector ki is always positive, all the negative 
terms are exactly canceled by positive terms. It can never been zero. □ 

Theorem 7 

It is impossible to have a collinear Four-Body central configuration with two 
charges of one sign on one side of the line and two charges of opposite sign 
on the other side. 

Proof: Assume charges 1 and 2 are positive and charges 3 and 4 are 
negative. The forces for this distribution of charges are 

(84) 
(85) 
(86) 
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The interior product of F with the vector k 2 is always positive definite. It 
can never been zero. □ 

Theorem 8 

It is impossible to have a collinear Four-Body central configuration with two 
charges of one sign inside the two charges of opposite sign. 
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Proof: Assume charges 1 and 4 are with different sign than charges 2 
and 3. The forces for this distribution of charges are 

F] mo vn-i rriA , . 

+ 7 , , , ( 88 ) 



miG a 2 (a + 6) 2 (a + 6 + c) 2 

F 2 mi m 3 m 4 

+ 177 ~ 7771—772" (89) 



m 2 G a 2 6 2 (6 + c) 2 
F 3 mi m 2 m 4 

777.3G (a + 6) 2 6 2 c 2 



(90) 



Fa mi mo m-? 

+ 7^7 i 7^ + -#- (91) 



m 4 G (a + 6 + c) 2 (6 + c) 2 c 

The interior product of F with the vector k3 is always positive definite, all 
the negative terms are exactly canceled by positive terms. It can never been 
zero. □ 

The possibility of collinear central configuration of gravitational charges 
of both signs is discovered by numerical computation of several cases, by 
means of the algorithm published in [H] using the same coordinate system. 
The possible particular permutation of coordinates is limited to a spherical 
triangle as represented in figure 2. Many central configuration with two par- 
ticles of one charge and the other with different charge have been computed 
when the sign of the charges alternates in the line: r\ and of one sign, r 2 
and ta of the opposite sign. Collinear central configurations with one charge 
different of the other has been computed in several cases when the different 
charge is located between two of the opposite charged particles. The pres- 
ence of one different charge perturb the positions of the three others which 
expand when the value of the different charge increases. Numerical examples 
will be e-mailed by request. 
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